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ABSTRACT 

We illustrate the description of correlated subsystems by studying the simple two-body Hydrogen atom. 
We study the entanglement of the electron and proton coordinates in the exact analytical solution. This 
entanglement, which we quantify in the framework of the density matrix formalism, describes correlations 
in the electron-proton motion. 
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1 Introduction 

Independent particle models often provide a good starting point to describe the physics of many-particle 
systems. In these models, the individual particles behave as independent particles that move in a potential 
field which accounts for the interaction with the other particles in an average sense. However, depending on 
the systems and the properties that arc studied, the results of the independent particle calculations arc not 
always sufficiently accurate. In that case it is necessary to correct for the fact that the motion of a single 
particle depends on the positions of the other particles, rather than on some average density. Consequently, 
in a system of interacting particles, the probability of finding two particles with given positions or momenta 
is not simply the product of the single particle probabilities: we say that the particles are 'correlated'. 

In a more general context, the problem can be formulated as the description of interacting subsystems (the 
single particle in the many-particle system being the example of the subsystem). A convenient theoretical 
framework to deal with this problem in quantum mechanics is provided by the density matrix theory, which 
is almost as old as quantum mechanics itself [|l|, The ever-present interest in density- matrix theory, in 
spite of its long history, is justified by the power of its description and the importance of its applications, 
which extends to the study of the very foundations of Quantum Theory ^, ||]. In this article, we briefly 
review some important aspects of density matrix theory and we illustrate its use for describing correlations 
of interacting subsystems by studying the simple, exactly solvable system of the Hydrogen atom. 

2 Subsystems and Quantum Correlations 

In very broad terms, giving the state of a physical system means providing the necessary information to 
evaluate all its observables quantitatively. The states of a quantum mechanical system are frequently rep- 
resented by vectors of unit norm in a Hilbert space Ti., in general of infinite dimension (hence their usual 
designation as 'state vectors'). If the system is in a state j^*), with (^'|^') = 1, the outcome of measurements 
of an observable quantity A has expectation value (A) equal to 

(A) = {^\A\^) , (1) 
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where we use a hat, 'A', to distinguish the operator representing the observable from the observable itself. 
If the vectors jwi), i — 1,2,... are a basis of the Hilbert space H, we can use the completeness relation, 
1 ~ X^m \''^m.){wm\i to Write the expectation value as 

(A) -^(^'|u;„0(i(;™|^|ii;„)(u;„|'J') . (2) 

m,n 

In matrix notation, the mn-matrix element of the A-operator in the |w)-basis is equal to Amn = {wm\A\wn) ■ 
The remaining ingredients of Eq. (H) can be collected as pnm = known as the density matrix 

10, 1^, ^, which is recognized as the matrix representing the so called density operator, p — in the 

chosen basis. The expectation value of Eq. (||) can then be written as 

(A) = y ^ A„inPnm j (3) 

which is the trace of the product of the A and p-matrices, (A) ~ Tr{Ap). In this expression, the state of 
the system is represented by the density matrix or, equivalently, by the density operator p. The form of this 
operator, together with normalization property of the underlying state vector j^*), lead to some important 
properties, which are of course shared by the corresponding density matrices. Namely, the density operator 
is a non-negative (meaning that {(f>\p\(f>) > for any vector self-adjoint operator (p = p^), of unit trace 
(Tr(p) = 1), and idempotent (meaning that p = p^, which follows from |\E')(\I'|\I')(*| = |^')(*|). 

The use of state vectors \'^) to describe the state of a quantum system is however not general enough to 
cover many frequently occuring situations. As we shall discuss in continuation, a suitable generalization is 
provided by the density matrix language, if only we relax the idempotency requirement. 

Often, when probing a system, only part of the total system is subjected to a measurement. For example, 
in a typical scattering experiment only the scattered particle is detected. The target system is left behind 
in a final state that could be different from its initial state, but this state is not observed directly. It is 
then natural to divide the total system into the target system and the single particle system of the scatterer. 
Also, in many experiments involving mesoscopic, or semi-macroscopic, quantum devices, it is very hard, if not 
impossible, to ensure efficient isolation from the environment, which then plays the role of a second subsystem 
coupled to the system of interest Q . Mathematically, such partitioning of the system into subsystems means 
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factorizing the Hilbert space 7i into the tensor product of the corresponding subspaces, 7i = Hu (S) If 
the Hu and Hy -spaces are spanned respectively by the bases \ui), i = 1, 2, ... and \vj), j = 1,2,..., then H 
is spanned by the product states \ui)\vj). Accordingly, a state vector \^) describing the state of the system 
can be expanded as 

m =Y,di,j\ui)\vj) . (4) 

Suppose now that one wishes to describe the state of one of the subsystems alone, say subsystem u. The 
observable quantities of this subsystem are represented by operators which act non-trivially on the vectors of 
Hu while acting on the vectors of Hy simply as the unit operator. The expectation value of such operators 
can be obtained from the reduced density matrix for the u-system, which is defined as the trace of the full 
density matrix p over the v-subspace. When the entire system is in the state described by l^*) this reads 

Pu = Try \^^\ui)\vj)dijdl^i{uk\{vi\ 

i,k j 

We can understand the trace as a sum J2j over the probabilities for the u-subsystem to be in any possible \vj)- 
state. For example, in the scattering experiment where only the state of the detected particle is determined, 
the final state of the target (i;)-system is not measured, and one has to sum over the probabilities of finding 
the target in all possible |i);)-states. 

An important limiting situation is that in which the overall state vector is itself a product of a 
u-vector and a w-vector, |^) = |u) \v). For these product states, the density matrix p = |^')(^'|, factorizes, 

p=\u)\v) {v\{u\=p^p- , (6) 

and the reduced density p^ is simply given by This is just of the form of a density operator associated 

with the state vector \u), which can then also be used to describe the state of the subsystem. In general, the 
state vector | \E') cannot be factorized in this fashion, however. In this case one says that the two subsystems 
are in an 'entangled' state. The density matrix of an entangled state does not factorize and appears as a 
series of u and w-products: the subsystems are correlated. The number of terms gives some indication of the 
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departure from the uncorrelated density matrix, but as the number of terms depends on the choice of basis 
in the u and w-spaces, this is not a good measure of correlation or entanglement. A way out of this difficulty 



identifies a 'natural' basis (in the sense of being determined by the structure of the entangled state itself) in 
which the description of the entanglement is achieved with maximum simplicity. 

The Schmidt basis is a product basis, in the usual sense that it is written as the tensor product of two 
particular bases, one for the Tiu space and another for the 7i„ space. In order to find these two bases one 
looks for the eigenvectors of the reduced density matrices. If the state is not entangled, |^) — \u)\v), and 

each have one single eigenvector of non-zero eigenvalue, \u) for p" and \v) for p", and the numerical value 
of the corresponding eigenvalues is equal to 1. One then completes the basis in each subspace by including 
enough additional orthonormal vectors which are moreover orthogonal to the single 'relevant' vector with 
nonzero eigenvalue. Each one of these aditional basis vectors is then an eigenvector of the corresponding 
reduced density with eigenvalue zero, and the large degree of arbitrariness in choosing them reflects the large 
degeneracy of the corresponding eigenvalue zero. If, on the other hand, the subsystems are correlated, then 

and p"" have a spectrum of nonvanishing eigenvalues. Interestingly, as we prove below, the eigenvalues 
of the and p^-operators are the same. To see that, we assume that the p"-matrix in the \ui) basis, is 
diagonalized by the unitary transformation U UiiU*j — Sij; J2i UiiU*^ — Sji), 



is however provided by a simple and remarkable result due to Schmidt [Q (see also |[ ^ ) which essentially 




(7) 
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The p"-matrix elements are given by Eq.(||), p"^. ~ J2j '^ij'^kj' ^iis^t Eq.(^ can be written as 




(8) 



Multiplying Eq.(^) by d*;/, and summing over i, we obtain 




(9) 



Defining Vj,i 



J2k d^jUki, we find with Eq.(^ that the vector of components {Vu, V21, ■ ■ ■ ,Vji, . . .) is an 
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eigenvector of with eigenvalue A;: 

Y^pl^V,i = \iViv , (10) 

where p^i^ = J^i'^u'^'ij- 1^ follows that and p" have the same eigenvalues A. If the nonvanishing eigen- 
values are not degenerate, then all these vectors are automatically orthogonal. And if, moreover, the set of 
eigenvectors with nonvanishing eigenvalue is still not complete, one still may complete the bases with ad- 
ditional orthogonal vectors which are again eigenvectors of the respective reduced densities with eigenvalue 
zero. Expanding j^*) in the product basis constructed in this way, which thus includes the eigenvectors \ux) 
of p" and I^a) of p^ ,1"^) = J2x c\\u\)\v\), the reduced density matrices take on the form 

= J2\cx\'\ux){u^\, (11) 

A 

which is diagonal and shows that the A-eigenvalues are equal to A = |cAp. This also shows that the 
eigenvectors with zero eigenvalue do not participate in the expansion. The eigenvalues A are thus both the 
probabilities of finding subsystem u in the states \u\) and the probabilities of finding subsystem v in the states 
\vx). The set of numbers |caP can therefore be interpreted as a distribution of occupation probabilities. In 
fact, the unit trace condition on the complete density operator p ensures that |caP = 1- For uncorrelated 
subsystems, there is just one nonvanishing probability and this distribution has a vanishing 'width' or 
standard deviation. It is then reasonable to use the standard deviation as a measure of the correlation 
pO| . Note that when the two subsystems are correlated, i.e., there is more than just one nonvanishing 
occupation probability |caP, the reduced density matrix Eq. ( [ll| ) is no longer idcmpotcnt, although all the 
other properties listed above for p still hold. As a result of this, the state of an entangled subsystem cannot 
be described in terms of a state vector in the corresponding Hilbert space, and one is forced to use the density 
matrix language. In this case one says that the subsystem is in a 'mixed' state. 'Pure' states, associated 
with definite state vectors, are on the other hand always associated with idempotent density operators. 

In summary, the foregoing discussion shows that a subsystem of a larger quantum system is, in general, 
in a mixed state, even if the state of the system as a whole is pure. Furthermore, the Schmidt analysis 
shows that in this case the state of the complementary subsystem has exactly the same degree of impurity, 
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in the sense that both reduced density matrices have the same spectrum. This corresponds to the impurity 
of both subsystems being due to their mutual entanglement in the given overall pure state. To the extent 
that a given quantum system is not really isolated, but interacts with other systems (possibly a non-descript 
'environment'), this analysis also shows that assuming the purity of its state is inadequately restrictive. In 
the following section we examine the spatially homogeneous, but possibly not pure, states of a single quantum 
particle, independently of any dynamical processes that may in fact give real existence to such states. An 
example of such a dynamical process will be given next. 

3 Single Particle System 

For the sake of illustration, we apply the density matrix description to a single, 'elementary' particle system. 
In spite of the extreme simplicity of this type of system, the density matrix language enables us to consider 
states that are not accessible within the standard state vector (or wavefunction) description - mixed states 
- which are of interest for the discussion of correlated many-particle systems. 

In the coordinate representation, which corresponds to the choice of the eigenfunctions |r) of the coordi- 
nate operator as basis, the density matrix p of a pure single particle system represented by the state vector 
1^), takes on the form 

p(r,r') = (r|*)(*|r') = *(r)**(r') , (12) 

where (r|^) = ^'(r) is the wavefunction. In the general case, the single particle density matrix is an object 
p{r, r') with the properties of hermiticity, p(r, r') = p*{r', r), unit trace, / d'^rp{r, r) = 1, and non-negativity, 
J d^r J (i^r'(/)*(r)p(r, r')(^(r') > 0. Note that idempotency is not being required in order to allow for mixed 
states. 

We shall be concerned with a particular class of states which we shall refer to as homogeneous states. 
By definition, these states are translationally invariant, meaning that translation by an arbitrary position 
vector R leaves the density matrix, and hence the state itself, unaltered: p{r + R, r' + R) = p{r, r'). As a 
consequence, all the density matrix elements depend only on relative positions, specifically, p(r, r') = p(r— r'). 
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In this case, the Fourier transform /5(k) of p{y — r'), 

P(r - r') = ^ j d?k p(k) exp(*k • [r - r']) , (13) 

plays a dual role in the description of the homogeneous state: (1) it is proportional to the momentum 
distribution of the system, and (2) it gives the eigenvalues of the density matrix. We shall now prove both 
assertions. 

Let us first evaluate the expectation value (p) of the momentum. If the system is pure and described by 
the quantum state \^), then 

(p) = j rf-V ^*{r) [-iW] l'(r) 

= -ih J (fr d^r' S{r -r')Vrp{r,r') . (14) 



The trace prescription for calculating (p), obtained above, is also valid for mixed states and inserting Eq.(13), 
we obtain 

(p> = -ih J d^r J dV(5(r - r')Vrp(r, r') 

= / "^^^ / / '^^^ [^^] '°(^) ^''P^*^ ■ " ^']) " (1^) 

= (Ijl / [^J^] ^(J^) ' 

where we have used a large quantization volume to avoid problems related to states of infinite norm. This 
normalization scheme, commonly called box normalization, assumes that the particle is confined to a large 
volume Q, which is taken to be infinite at the end of the calculation, — s- cx). A plane wave state is then 
normalized as exp(ik • r)/y^. Substituting p = Tik in the integration of Eq.(p^, we find that 



d^p p /(p) (16) 

where /(p) = [fl/{2TTh)^]p{p/h). Similarly, the expectation value of higher-order powers of momentum 
components are equal to the corresponding higher-order moments of the /(p)-function, e.g. ((p • x)^) — 
J d^p (p • x)^ /(p)- Thus, /(p) may be interpreted as a momentum distribution function. 
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To see that p(k) gives the eigenvalues of the density matrix, we use Eq. (^3|) to evaluate 

which we recognize as the eigenvalue equation of the density matrix (Eq.(^) in coordinate representation. 
From Eq. , it follows that the eigenvectors of the density matrix are the plane wave states of momentum 
k and that the corresponding eigenvalues are equal to p(k). 

The density matrix of a homogeneous single particle system does generally not describe a pure system. 
To see that, we remind the reader that a pure system is characterized by an idempotent density matrix, 

= p. For the homogeneous single particle system, this implies that the eigenvalues of the density matrix 
satisfy p(k)^ = /5(k). Consequently, for a pure homogeneous single particle system, we find that p(k) = or 
p(k) = 1. Furthermore, the unit trace condition on the density matrix implies 

where in the last step we reconverted the momentum integral to the sum of discrete momenta appropriate 
to the adopted volume quantization. Thus, the homogeneous single particle system that is pure can only 
be described by a density matrix with a single eigenvalue, say p(k.i), equal to 1, and all other eigenvalues 
p(k) = 0,k 7^ ki. Such density matrix represents a particle with a definite momentum (zero standard 
deviation of the distribution of density matrix eigenvalues), corresponding to a plane wave wavefunction, 
(rj^*) = exp(iki ■r)/\^. In contrast, an impure single particle system, has a spectrum p(k) of density matrix 
eigenvalues. The standard deviation of the p(k)-distribution is a measure of the correlation or entanglement 
of the single-particle system . 

It should be remarked, finally, that all the above discussion had a purely kinematical character, in the 
sense that we dealt with possible states of a quantum mechanical system and with their inherent obervable 
properties, without regard to any dynamical processes that might have led to their preparation. In the 
following section we consider a simple but definite composite dynamical system and use the tools developed 
here in order to study the entanglement effects generated by definite interaction processes. 
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4 Hydrogen Atom 

We now apply the density matrix formahsm to stationary states of a system that can display correlation: the 
two-body hydrogen atom. The hydrogen atom consists of an electron and a proton, interacting by means of 
the Coulomb potential. The hamiltonian of the electron-proton system is 

H=-2^ + -^--^ , (19) 

2mo 2mp |rc-rp| 

where p and r denote the momentum and position operators in coordinate space, and the e and p-subscripts 
indicate the electron and the proton. 

In fact, the hydrogen atom can be treated as two uncorrelated subsystems. Transforming to relative, 
r = Tp—Ve, and center-of-mass coordinates, R = {meVe+mprp) / AI , where M is the total mass, M = me+mp, 
the hamiltonian separates into a center-of-mass term which takes on the form of a free-particle hamiltonian, 
^^cm = Pfl/SM, and an 'internal' hamiltonian which governs the relative motion of the electron and the 
proton, Hiat — p'l/'^rar — e^/r, where rrir is the reduced mass, — m^^ + rrip^. The product of Hem and 
-ffint-eigenstates, 

*(R,r) = 0int(r)exp(i^-R)/Vo (20) 

is an eigenstate. According to the definitions introduced in section 2, the internal and center-of-mass sub- 
systems for the state ( po|) are not correlated. If the hydrogen atom is in its atomic ground state, the internal 
wavefunction, 0int, is equal to 

(^i„t(r) = — exp(-r/ao) (hydrogen m Is - state) , (21) 

VTT \ao/ 

where ag is the Bohr-radius, gq = y. 



Partitioning the hydrogen atom differently into electron and proton subsystems gives an entangled state: 
nre,rp) = ^Ure - r^) exp • Kre + m^r^] ^ ^^2) 
The reduced density matrix for the electron system can be obtained by taking the trace over the proton 
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basis set of coordinate eigenfunctions: 



meP 

hM 



^^^P ( ■ [""e ~ I"e] ) Pint(l"e, Tg) (23) 



where we introduce the electron density matrix pi^t for electrons in atoms at rest (P = 0). It is interesting 
to note that the electron subsystem of the hydrogen atom is homogeneous. To prove that, it is sufficient to 
show that Pint depends on y^. and r'g as — v'^. Substituting rj. — rp = y and writing the argument of (/)i„t , 
Te — Tp, as Te — r'g + y, we do indeed find that 

d^rp 0i„t (r - rp)0*„t (r^ - rp) = j d^y 0int (re - r', + y)(f>l, (y) , (24) 



which depends solely on — r'g. Consequently, the reduced electron density matrix (^3|) only depends on 
the relative position, indicating a homogeneous subsystem. This might appear surprising: we do not think 
usually of the electron in hydrogen as a homogeneous system. However, the reduced density matrix describes 
the observation of hydrogenic electrons in the assumption that the proton is 'invisible'. Furthermore Eq.(|22 
describes a hydrogen of fixed momentum (e.g. a beam of hydrogen atoms of well-defined velocity) so that 
the center-of-mass position of the atom is undetermined. It is then equally likely to observe the electron in 
any position, and the electron subsystem is homogeneous. 

Consequently, as discussed in the previous section, the Fourier-transform of the density-matrix plays a 
central role. The density-matrix pint of a Is-electron in a hydrogen atom at rest (P = in Eq.(p3|)), has a 
simple analytical Fourier-transform, 

1 GAnal 

The Fourier transform p for an electron in an atom of arbitrary momentum P is a translation in k-space of 
Pint- To show that, we insert the inverse Fourier transform of p in the expression for the density matrix (|2^), 

P(re,0 = ^ /rf'fc'A„t(k')exp(z[^+k'] .[r,-r^]) . (26) 
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The substitution, k = k' + ^^f^, then leads to the Fourier-transform: 



p(r„ r^) = ^ Jd'k (k - ^) exp(zk • [r, - r^]) , (27) 
from which it follows that 

p(k)=pi„t(k-^) . (28) 

Thus, compared to pi^t, p is displaced in k-space by . Similarly, the momentum distribution /(p) = 
[fl/{2iTh)^]pint{p/K), is related to the electron momentum distribution /int(p) for atoms at rest, by means 
of a simple displacement in momentum space: 

/(p)=/int(p-^P) , (29) 

where 

(27ra)3l + (£f) ■ ^^^> 
Equation ( p9| ) expresses a Gallilean transformation from the center-of-mass reference frame (P — 0) to the 
lab frame of reference in which the atom is moving with velocity P/M. If the electron in the center-of- 
mass frame has velocity v — p/rrie, then, in the lab frame, its velocity is observed as v' = p/rn-e — P/M, 
corresponding to a momentum equal to mgv' = p — -^P- 

Since / and /int are equal up to a translation in momentum space, their widths, or standard deviations, 



are equal. In the center-of-mass frame, (p) = 0, so that the standard deviation Ap = \J (p^) — ((p))^ is 
equal to 



Ap=Jj (Ppp^fUv) 



= - . (31) 
ao 

Thus, although the electron subsystem is homogeneous, as a consequence of the electron-proton correlation, a 
measurement of the electron momentum can yield a finite range of values. In accordance with the Heisenberg 
uncertainty principle (the correlations confine the electron in the center-of-mass frame to a region of size 
~ ao) the size of this region in momentum space, or more precisely, the standard deviation Ap of the 
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momentum distribution, is equal In the context of describing correlating systems, Ap is a measure of 
the strength of the correlation: the stronger the electron and proton are correlated, the higher the value of 
ao and the larger the region in momentum space (Ap) over which the electron momentum is 'spread out'. 

5 Conclusion 

In this paper, we have reviewed aspects of density matrix theory which pertain to the description of the 
entanglement of correlated subsystems. The mutual entanglement of subsystems of a larger (pure) system is 
conveniently quantified by the standard deviation of the reduced density matrix eigenvalues of the entangled 
subsystems. For the simple (but relevant) case of homogeneous single particle subsystems, the density matrix 
eigenvalues give the particle's momentum distribution. Thus, the standard deviation of the momentum 
distribution is a measure of the entanglement of the single particle system, as we have illustrated for the 
specific example of the electron in the two-body hydrogen atom. 
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